Fourth order gravity and experimental constraints on Eddington parameters 
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PPN-limit of alternative theories of gravity represents a still controversial matter of debate and 
no definitive answer has been provided, up to now, about this issue. By using the definition of the 
PPN-parameters 7 and /3 in term of f(R) theories of gravity, we show that a family of third-order 
polynomial theories, in the Ricci scalar R, turns out to be compatible with the PPN-limit and the 
deviation from General Relativity, theoretically predicted, can agree with experimental data. 

PACS numbers: 04.50.+h, 04.25.Nx, 98.80.-k 

o : 

Keywords: Extended theories of gravity; PPN limit; experimental tests of general relativity. 

1. General Relativity (GR) is the cornerstone theory among the several attempts proposed to describe gravity. It 
represents an elegant approach furnishing several phenomenological predictions and its validity, in the Newtonian 
limit regime, is experimentally probed up to the Solar System scales. However, also at these scales, some conundrums 
come out as the indications of an apparent, anomalous, long-range acceleration revealed from the data analysis of 
I Pioneer 10/11, Galileo, and Ulysses spacecrafts which are difficult to be framed in the standard scheme of GR and its 
low energy limit 0,0- Furthermore, at galactic distances, huge bulks of dark matter are needed to provide realistic 
t-H ' models matching with observations. In this case, retaining GR and its low energy limit implies the introduction of an 
actually unknown ingredient. We face a similar situation even at larger scales: clusters of galaxies are gravitationally 
stable and bound only if large amounts of dark matter are supposed in their potential wells. Finally, an unknown 
form of dark energy is required to explain the observed accelerated expansion of cosmic fluid. Summarizing, almost 
95% of matter-energy content of the universe is unknown in the framework of Standard Cosmological Model while 
we can experimentally probe only gravity and ordinary (baryonic and radiation) matter. Considering another point 
— «. . of view, anomalous acceleration (Solar System), dark matter (galaxies and galaxy clusters), dark energy (cosmology) 
could be nothing else but the indications that shortcomings are present in GR and gravity is an interaction depending 
1 . on the scale. The assumption of a linear Lagrangian density in the Ricci scalar R for the Hilbert-Einstein action could 
^Jy be too simple to describe gravity at any scale and more general approaches should be pursued to match observations. 
. . , Among these schemes, several motivations suggest to generalize GR by considering gravitational actions where generic 
^ • functions of curvature invariants are present. Specifically, actions of the form 



^ ! A = I d A x y 



f(R) + £r, 



(1) 



where f(R) is an analytic function of R and C m is the standard matter Lagrangian density, result particularly 
interesting. The variation of Q gives rise to fourth-order field equations 

G aP = R a p - \g aP R = T$T + T^/f'(R) (2) 
where the curvature stress - energy tensor is defined as 

T^ rv = i~g a p [f(R) - Rf'(R)] + f(R)^(g a ^ - 5^.9^)} > (3) 

with prime denoting derivative with respect to R and T™p the standard matter contribution. For f{R) — R, the 
curvature stress - energy tensor identically vanishes and Eqs.© reduce to the standard second -order Einstein field 
equations. From Eq.Q, it is clear that the curvature stress - energy tensor plays the role of a further source term in 
the field equations so that it can be considered as an effective fluid of purely geometric origin. 

This approach is physically motivated by several unification theories of fundamental interactions and by the field 
quantization on curved space-times Q. At cosmological level, it is well known that further curvature contributions 
can solve shortcomings at early epochs (giving rise to inflationary solutions and explaining the today observed 
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accelerated behavior by a sort of curvature quintessence, i.e. curvature can act as a fluid implementing acceleration 
II 110, HUGH This result can be achieved in metric and affine (Palatini) approaches mElEQIllHIllE 
[jjj, |2£j ■ In addition, reversing the problem, one can reconstruct the form of the gravity Lagrangian by observational 
data of cosmological relevance through a "back scattering" procedure. All these facts suggest that the function f(R) 
should be more general than the linear Hilbert-Einstein one implying that higher order gravity could be a suitable 
approach to solve GR shortcomings without introducing mysterious ingredients as dark energy and dark matter (see 

e.g. pm. 

In recent papers, some authors have confronted this kind of theories even with the Post Parameterized Newtonian 
prescriptions in metric and Palatini approaches. However, the results seem controversial since in some cases [2^, 01 
it is argued that GR is always valid at Solar System scales and there is no room for other theories; nevertheless, 
some other studies (25[ [2(J find that recent experiments as Cassini and Lunar Laser Ranging allow the possibility 
that extended theories of gravity could be seriously taken into account. In particular, it is possible to define PPN- 
parameters in term of f(R) functions and several classes of fourth order theories result compatible with experiments 
in Solar System pfij . 

In this letter, we follow a different approach. Starting from the definitions of PPN-parameters in term of a generic 
analytic function f(R) and its derivatives, we deduce a class of fourth order theories, compatible with data, by means 
of an inverse procedure which allows to compare PPN-conditions with data. As a matter of fact, it is possible to show 
that a third order polynomial, in the Ricci scalar, is compatible with observational constraints on PPN-parameters. 
The degree of deviation from GR depends on the experimental estimate of PPN-parameters. 

A useful method to take into account deviation with respect to GR is to develop expansions about the GR solutions 
up to some perturbation orders. A standard approach is the Parameterized-Post-Newtonian (PPN) expansion of the 
Schwarzschild metric. In isotropic coordinates, it is 

ds2 = (l^f ) a<ft2 - (l + j)V 2 +r~W) (4) 

where r g — 2GM/c 2 is the Schwarzschild radius. Eddington parameterized deviations with respect to GR, considering 
a Taylor series in term of r g /r assuming that in Solar System, the limit r g /f) <C 1 holds [2?|. The resulting metric is 



ds 2 
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dt 2_i 1 + 7 !l + ...)Uf 2 +fW) (5) 



where a, f3 and 7 are unknown dimensionless parameters (Eddington parameters) which parameterize deviations with 
respect to GR. The reason to carry out this expansion up to the order {r g /r) 2 in 500 and only to the order (r g /f) 
in gij is that, in applications to celestial mechanics, g^ always appears multiplied by an extra factor v 2 ^ (M/f). It 
is evident that the standard GR solution for a spherically symmetric gravitational system in vacuum is obtained for 
a = (3 = 7 = 1 giving again the Schwarzschild solution. Actually, the parameter a can be settled to the unity due to 
the mass definition of the system itself jS^. As a consequence, the expanded metric JSJ) can be recast in the form : 



ds 2 
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dr 2 - r 2 dil 2 . 
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where we have restored the standard spherical coordinates by means of the transformation r — r (l + ?§) . The 
two parameters /3, 7 have a physical interpretation. The parameter 7 measures the amount of curvature of space 
generated by a body of mass M at radius r. In fact, the spatial components of the Riemann curvature tensor are, at 
post-Newtonian order, 

Rijki = %r£N m (7) 



independently of the gauge choice, where Nijki represents the geometric tensor properties (e.g. symmetries of the 
Riemann tensor and so on). On the other side, the parameter f3 measures the amount of non-linearity (~ (r g /r) 2 ) in 
the goo component of the metric. However, this statement is valid only in the standard post-Newtonian gauge. 

If one takes into account a more general theory of gravity, the calculation of the PPN-limit can be performed 
following a well defined pipeline which straightforwardly generalizes the standard GR case [23. A significant devel- 
opment in this sense has been pursued by Damour and Esposito-Farese [2^, El which have approached to the 
calculation of the PPN-limit of scalar-tensor gravity by means of a conformal transformation g^ v = F{<j))g^ v to the 
standard Einstein frame. In fact, a general scalar-tensor theory 



A = 



F{4>)R+\<t>;^ -V{<l>) + L r , 



(8) 
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TABLE I: Summary of the three approaches: Scalar- Tensor (ST), Einstein +ip, (E + <f>) and f{R) and their relations at 
Lagrangians, field equations and solutions levels. The solutions are in the Einstein frame for the minimally coupled case while 
they are in Jordan frame for f(R) and ST-gravity. Clearly, f(R) and ST theories can be rigorously compared only recasting 
them in the Einstein frame. 



where V(4>) is the self-interaction potential and F(4>) the non-minimal coupling, can be recast as : 

A = I 
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where 



dtp 



d\nF(<j>) 



2F{4>) ' 



(9) 
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8ttG = 1, V(ip) = F 2 ((j>)V(<fi), C m = F 2 (cj))£ m . The first consequence of such a transformation is that now the 
non-minimal coupling is transferred to the ordinary-matter sector. In fact, the Lagrangian C m is dependent not only 
on the conformally transformed metric and matter fields but it is even characterized by the coupling function 
F(<f)) 32] . This scheme provides several interesting results up to obtain an intrinsic definition of 7, (3 in term of 
the non-minimal coupling function F((j>). The analogy between scalar-tensor gravity and higher order theories of 
gravity has been widely demonstrated [33L l34l l35l |. Here an important remark is in order. The analogy between 
scalar-tensor gravity and fourth order gravity, although mathematically straightforward, requires a careful physical 
analysis. Recasting fourth-order gravity as a scalar-tensor theory, often the following steps, in terms of a generic 
scalar field ip, are considered 



f(R) + £r 



F'(iP)R + F(4>) - F'(V# + Cm -> F'{iP)R - V(i>) + £, 



(11) 



where, by analogy, ip — > R and the "potential" is V(ip) — F(tp) — F'(ip)ip. Clearly the kinetic term is not present 
so that (fTTfl is usually referred as a Brans-Dicke description of f(R) gravity where (jJbd = 0. This is the so-called 
O'Hanlon Lagrangian |3^|. However, the typical Brans-Dicke action is 



A = 



IpR - UJBD 



1> 
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where no scalar field potential is present and lobd is a constant. In summary, O'Hanlon Lagrangian has a potential but 
has no kinetic term, while Brans-Dicke Lagrangian has a kinetic term without potential. The most general situation 
is in JSJl where we have non-minimal coupling, kinetic term, and scalar field potential. This means that fourth-order 
gravity and scalar tensor gravity can be "compared" only by means of conformal transformations where kinetic and 
potential terms are preserved. In particular, it is misleading to state that PPN-limit of fourth order gravity is bad 
defined since these models provide ojbd — and this is in contrast with observations |23ll24j . 

Scalar-tensor theories and f(R) theories can be rigorously compared, after conformal transformations, in the Ein- 
stein frame where both kinetic and potential terms are present. With this consideration in mind, F(<f>) and f'(R) can 
be considered analogous quantities in Jordan frame and then the PPN limit can be developed 1 . 

Starting from this analogy, the PPN results for scalar-tensor gravity can be extended to fourth order gravity |25l. 
In fact, identifying (f> — * R 35], it is possible to extend the definition of the scalar-tensor PPN-parameters |2Sl.l36| to 
the case of fourth order gravity : 



1 = -- 



f"(RY 



f'(R) + 2f"(RY 



0-1 = 7 



f'(R)-f"(R) \ ch 
2f'(R) + 3f"(R) 2 dR- 



(13) 



To be precise, conformal transformations should be operated "before" performing PPN-limit and results discussed in the same frame. 
A back conformal transformation, after PPN limit, could be misleading due to gauge trobles. 
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In [25j , these definitions have been confronted with the observational upper limits on 7 and (3 coming from Mercury 
Perihelion Shift [32J and Very Long Baseline Interferometry (38|. Actually, it is possible to show that data and 
theoretical predictions from Eqs. 1)13(1 agree in the limits of experimental measures for several classes of fourth order 
theories. Such a result tells us that extended theories of gravity are not ruled out from Solar System experiments 
but a more careful analysis of theories against experimental limits has to be performed. A possible procedure could 
be to link the analytic form of a generic fourth order theory with experimental data. In fact, the matching between 
data and theoretical predictions, found in |25j . holds provided some restrictions for the model parameters but gives 
no general constraints on the theory. In general, the function f(R) could contain an infinite number of parameters 
(i.e. it can be conceived as an infinite power series [34]]) while, on the contrary, the number of useful relations is finite 
(in our case we have only two relations). An attempt to deduce the form of the gravity Lagrangian can be to consider 
the relations 113fl as differential equations for f(R), so that, taking into account the experimental results, one could 
constrain, in principle, the model parameters by the measured values of 7 and (3. 

3. The idea is supposing the relations for 7 and (3 as differential equations. This hypothesis is reasonable if the 
derivatives of f(R) function are smoothly evolving with the Ricci scalar. Formally, one can consider the r.h.s. of the 
definitions (|13|l as differential relations which have to be matched with values of PPN-parameters. In other words, 
one has to solve the equations (|13fl where 7 and (3 are two parameters. Based on such an assumption, on can try to 
derive the largest class of f(R) theories compatible with experimental data. In fact, by the integration of Eqs. i|13[l . 
one obtains a solution parameterized by (3 and 7 which have to be confronted with the experimental quantities (3 exp 
and j exp . 



Assuming f'(R) + 2/"(i?) 2 ^ and defining A 
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we obtain from l|13|) a differential equation for f(R): 



27-I 
[f"(R)} 2 = Af'(R) . 

The general solution of such an equation is a third order polynomial f(R) = aR 3 + bR 2 + cR H 

have to satisfy the conditions : a = b = c = and d ^ (trivial solution) or a = yj, b = ±^-^, 
the general solution for the non-trivial case, in natural units, reads 



(14) 

d whose coefficients 
with c, d =/= 0. Thus, 
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It is evident that the integration constants c and d have to be compatible with GR prescriptions and, eventually, 
with the presence of a cosmological constant. Indeed, when 7 — > 1, which implies f(R) — > cR + d, the GR- limit is 
recovered. As a consequence the values of these constants remain fixed (c=l and d = A, where A is the cosmological 
constant). Therefore, the fourth order theory provided by Ea. (|15f) becomes 
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where we have formally displayed the two branch form of the solution depending on the sign of the coefficient entering 
the second order term. Since the constants a, b, c, d of the general solution satisfy the relation 3ac — b 2 = 0, one can 
easily verify that it gives : 



dj 
dR 



f»(RY 



f±(R) 



dRf>(R)+2f"(Rf 



f±(R) 



(17) 



where the subscript / ± (_r) refers the calculation to the solution l|ltj|) . This result, compared with the second differential 
equation Eq. (|13l) . implies 4(/3 — 1) = , which means the compatibility of the solution even with this second relation. 
4. Up to now we have discussed a family of fourth order theories (|16H parameterized by the PPN-quantity 7; on the 
other hand, for this class of Lagrangians, the parameter j3 is compatible with GR value being unity. 

Now, the further step directly characterizes such a class of theories by means of the experimental estimates of 7. In 
particular, by fixing 7 to its observational estimate "fexp, we will obtain the weight of the coefficients relative to each 
of the non-linear terms in the Ricci scalar of the Lagrangian (|16fl . In such a way, since GR predictions require exactly 
"fexp — Pexp = 1 , in the case of fourth order gravity, one could to take into account small deviations from this values 
as inferred from experiments. Some plots can contribute to the discussion of this argument. In Fig^ the Lagrangian 
(|16ll is plotted. It is parameterized for several values of 7 compatible with the experimental bounds coming from 
the Mercury perihelion shift (see Table. 1 and The function is plotted in the range R > 0. Since the property 

f+(R) = —f-(—R) holds for the function (|16|) . one can easily recover the shape of the plot in the negative region. As 
it is reasonable, the deviation from GR becomes remarkable when scalar curvature is large. 
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Mercury Perihelion Shift 


\2j-/3- 1| < 3 x 1(T 3 


Lunar Laser Ranging 


4/3-7-3 = -(0.7 ±1) x 1(T 3 


Very Long Baseline Interf. 


|7 - 1| = 4 x 10" 4 


Cassini Spacecraft 


7-1 = (2.1 ±2.3) x 10~ 5 



TABLE II: A schematic resume of recent experimental constraints on the PPN-parameters. They are the perihelion shift of 
Mercury |3T|. the Lunar Laser Ranging 13(1 . the upper limit coming from the Very Long Baseline Interferometry [3H and the 
results obtained by the estimate of the Cassini spacecraft delay into the radio waves transmission near the Solar conjunction 
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FIG. 1: Plot of the two branch solution provided in Ea. l)lrift . The f+(R) (dotted line) branch family is up to GR solution 
(straight line), while the one indicated with f-(R) (dotted-dashed line) remains below this line. The different plots for each 
family refer to different values of 7 fulfilling the condition 7 — 1| < 10 -4 and increased by step of 10 -5 . 

In order to display the differences between the theory (|16|l and Hilbert-Einstein one, the ratio f(R)/R is plotted in 
Fig|21 Again it is evident that the two Lagrangians differ significantly for great values of the curvature scalar. It is 
worth noting that the formal difference between the PPN-inspired Lagrangian and the GR expression can be related 
to the physical meaning of the parameter 7 which is the deviation from the Schwarzschild-like solution. It measures 
the spatial curvature of the region which one is investigating, then the deviation from the local flatness can be due 
to the influence of higher order contributions in Ricci scalar. On the other hand, one can reverse the argument and 
notice that if such a deviation is measured, it can be recast in the framework of fourth order gravity, and in particular 
its "amount" indicates the deviation from GR. Furthermore, it is worth considering that, in the expression (|16|l . the 
modulus of the coefficients in 7 (i.e. the strength of the term) decreases by increasing the degree of R. In particular, 
the highest values of cubic and squared terms in R are, respectively, of order 10 -4 and 1CP 2 (see Fig[3J) then GR 
remains a viable theory at short distances (i.e. Solar System) and low curvature regimes. 
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FIG. 2: The ratio f(R)/R. It is shown the deviation of the fourth order gravity from GR considering the PPN-limit. Dotted 
and dotted-dashed lines refer to the f+(R) and f-(R) branches plotted with respect to several values of 7 (the step in this case 
is 2.5 x 10~ 5 ). 
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FIG. 3: Plot of the modulus of coefficients: a (R s ) (line) and b (R 2 ) (dashed line). The choice of plotting the modulus of 
coefficients is the consequence of two solution for f(R)- 

A remark is in order at this point. The class of theories which we have discussed is a third order function of the 
Ricci scalar R parameterized by the experimental values of the PPN parameter 7. In principle, any analytic f(R) 
can be compared with the Lagrangian provided suitable values of the coefficients. However, more general results 
can be achieved relaxing the condition [3=1 which is an intrinsic feature for (|16f) (see for example |25j). These 
considerations suggest to take into account, as physical theories, functions of the Ricci scalar which slightly deviates 
from GR, i.e. f(R) = JqR( 1+ ^ with e a small parameter which indicates how much the theory deviates from GR 
|4l| . In fact, supposing e sufficiently small, it is possible to approximate this expression 



This relation can be easily confronted with the solution (|16fl since, also in this case, the corrections have very small 
"strength" 0. 

5. We have shown how a polynomial Lagrangian in the Ricci scalar R, compatible with the PPN-limit, can be 
recovered in the framework of fourth order gravity. The approach is based on the formulation of the PPN-limit of 
such gravity models developed in analogy with scalar-tensor gravity 25] . In particular, considering the local relations 
defining the PPN fourth order parameters as differential expressions, one obtains a third-order polynomial in the Ricci 
scalar which is parameterized by the PPN-quantity 7 and compatible with the limit (3=1. The order of deviation 
from the linearity in R is induced by the deviations of 7 from the GR expectation value 7 = 1. Actually, the PPN 
parameter 7 may represent the key parameter to discriminate among relativistic theories of gravity. In particular, this 
quantity should be significatively tested at Solar System scales by forthcoming experiments like LATOR |43j. From 
a physical point of view, any analytic function of R, by means of its Taylor expansion, can be compared with l|16fl . 
Therefore, a theory like f(R) = foR^ 1+<i \ indicating small deviations from standard GR, is in agreement with the 
proposed approach, so, in principle, the experimental 7 could indicate the value of the parameter e. In conclusion, 
one can reasonably state that generic fourth-order gravity models could be viable candidate theories even in the 
PPN-limit. In other words, due to the presented results, they cannot be a priori excluded at Solar System scales. 
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